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ABSTRACT 

General relativity predicts that two freely counter-revolving test particles in the exterior 
field of a central rotating mass take different periods of time to complete the same full orbit; 
this time difference leads to the gravitomagnetic clock effect. The effect has been derived for 
circular equatorial orbits; moreover, it has been extended via azimuthal closure to spherical 
orbits around a slowly rotating mass. In this work, a general formula is derived for the main 
gravitomagnetic clock effect in the case of slow motion along an arbitrary elliptical orbit in 
the exterior field of a slowly rotating mass. Some of the implications of this result are briefly 
discussed. 

PACS: 04.20. Cv; 04.80 
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1. Introduction 

The gravitomagnetic clock effect is basically a reflection of the fact that according to gen- 
eral relativity there is a special temporal structure around a rotating mass. More specifically, 
let us consider circular geodesies in the equatorial plane of a Kerr black hole of mass M and 
angular momentum J. Let t+ be the period of prograde (retrograde) motion along such 
an orbit according to asymptotically static inertial observers; then, t+ — t- = 47r J/(Mc^). 
Moreover, let r± represent the corresponding proper periods according to comoving clocks; 
then r+ — r_ ~ 47r J/(Mc^) for orbits with radius r >> 2GM/c^. It is remarkable that at 
0{c~^) this result is independent of the gravitational constant G and the radius of the orbit 
[1,2]. Various theoretical aspects of this effect have been investigated [3-10]. On the obser- 
vational side, the possibility of its detection has been considered by a number of authors 
[11-17]. 

The general relativistic calculation of the gravitomagnetic clock effect for a general orbit 
is quite complicated. Other than equatorial circular orbits, only spherical orbits have been 
considered thus far to first order in the rotational perturbation within the post-Schwarzschild 
approximation scheme [3] . To ameliorate this situation, we recently showed that it is possible 
to recover the main general relativistic results at the lowest order, i.e. 0(c~^), by considering 
the linear approximation of general relativity for the exterior geometry of a slowly rotating 
mass [18]. The equations of motion of a test particle in such a stationary field reduce, in the 
slow-motion approximation, to an equation of the Lorentz form, dv/dt — 'Eg+£g+{v/c)x'Bg, 
where E^, = —GMv/r^ is the Newtonian gravitoelectric field, £g is the post-Newtonian 
gravitoelectric field and B^,, 
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Jr^ - 3(J • r)r , (1) 



is the post-Newtonian gravitomagnetic field (see [19, 20] for reviews of gravitoelectromag- 
netism). It should be noted that there is a certain arbitrariness in the definition of the 
gravitoelectromagnetic field depending on how the analogy with electrodynamics is devel- 
oped; in particular, the present treatment is different from the development in [20]. In this 
connection, we remark in passing that the approach that is closest to electromagnetism is 
based on the gravitational Larmor theorem (cf. [3,7,20] and references therein). The explicit 
form of £g can be simply derived using the isotropic form of the Schwarzschild metric along 
the lines of equation (10) of [20]. To 0{c''^), Eg is given by c^Sy = (u^ - 4$) - A{Ey ■ v)v, 
where $ = GM/r is the Newtonian potential. Note that at a given point r along the orbit, 
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reversing the direction of motion, v — > — v, leaves Sg invariant while the gravitomagnetic 
perturbation changes sign. The contribution of £g to the orbital period is thus independent 
of the sense of orbital revolution. In fact, the orbital perturbation due to £g is planar and 
leads to the Einstein pericenter precession. Calculating the period of motion in this case 
via azimuthal closure (i.e. the time that it takes for the particle to go from = 0o to 
= 00 + 27r or vice versa), one finds that there is no clock effect when otherwise identical 
clockwise and counterclockwise orbital periods are compared. Hence, in the calculation of 
the clock effect to 0(c~^), the post-Newtonian gravitoelectric perturbation due to the mass 
M may be neglected. Therefore, for the sake of simplicity we essentially ignore £g in the 
following sections but return to it in the Appendix, where the gravitoelectric contribution 
to the orbital period is explicitly given. 

The gravitomagnetic clock effect has not yet been detected observationally; therefore, 
0(c~^) terms are the only ones that are of current interest. At this level of approximation, 
we derive the general effect in this paper for an eccentric orbit. Though our approach is 
general, we have in mind for specific applications the motion of artificial satellites around 
the Earth. In the absence of the gravitomagnetic force, the orbit would ideally be a Kep- 
lerian ellipse; therefore, we treat the gravitomagnetic force as a linear perturbation of the 
dominant Newtonian gravitoelectric force. Since such a perturbation cannot be "turned off" 
for motion around the Earth, the elements of the unperturbed elliptical orbit play the role 
of subsidiary variables in our treatment: They are not all directly measurable. The unper- 
turbed orbital plane, however, is given by the instantaneous plane of the osculating ellipse 
when the observations begin by definition at t = 0. 

In Section 2 we derive the general form of the perturbed orbit within our approximation 
scheme. The "period" of the motion is calculated in Section 3 using the notion of azimuthal 
closure. Section 4 presents the general formula for the gravitomagnetic clock effect to 0{c~^) 
as well as a discussion of some of its main implications. 



2. Eccentric orbit 

Let us consider a solution of the equations of motion 



GMr 2G 

T- + 



3(r • J)r X V + r v x J 



(2) 



which would correspond, in the absence of the perturbation, to the unperturbed elliptical 
orbit as in Figure 1. To simplify the analysis, we introduce the "inertial" coordinate system 
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(X, Y, Z) adapted to the unperturbed orbit (as in Fig. 1) such that the unperturbed orbit 
is given by X — po cos (fi,Y — po sin (f and Z — 0, where 

^ — ^ ^ (-^ _ eg sin -a) - - Co sin 'a)t=o ■ (3) 

1 + eo cos V 

Here Oq, eo,v = ip — go and m are respectively the semimajor axis, eccentricity, true anomaly 
and the eccentric anomaly of the unperturbed orbit; moreover, the angle go is the argument 
of pericenter as depicted in Figure 1. The Keplerian frequency and the period of this orbit are 
given by loq — {GM / a^Y^"^ and Tq = 2t^/ujq, respectively. Let us note for future reference that 
Po = ao(l — eocosit) and ao(l — eQ)po = eoi^osin'O, where an overdot represents differentiation 
with respect to time t and Lq — ^jGMaQ{l — Cq) is the specific angular momentum of the 
unperturbed orbit. 

To write the equations of motion in a convenient form, we express equation (2) in terms 
of the (X, y, Z) coordinate system. It proves useful to introduce the cylindrical coordinate 
system (p, (/?, Z) as well, where X = p cos (/? and Y = p sin if. In the "inertial" (X, Y, Z) 
system, the equations of motion can be expressed as X + GMX/ p^ = Fx and similarly 
for the Y and Z components. The advantage of this system is that along the orbit Z can 
be treated only to first order in the perturbation; therefore, the radial position of the test 
particle is in effect given by p = (X^ + y^)-*^/^ and the perturbing acceleration {Fx, Fy, Fz) is 
simply evaluated along the unperturbed orbit. The equations of motion can then be written 
in the (p, </?, Z) system as 

p-pv + —r-^p ' (4) 

p(p + 2p0^F^ , (5) 

Z+^Z = Fz . (6) 
P 

Here Fp = Fx cos (p + Fy sin (p and F^ = —Fx sin ip + Fy cos (p. The perturbing accelerations 
are given by Fp = eLo/p^, F^ = -epo/Po and 

^ 2LoGJsini .^.^ ^. . . . i 

^z^^r^ — 7\ 2T [2(1 + eo cost") sm(p + Co smf cos (pj , (7) 

where e = 2GJ cos and v — (p — g^. 

To solve the equations of motion to first order in e, we start with equation (5) and write 
it in the form d{p'^(p)/dt — pF^, which expresses the rate of change of the Z-component 
of the orbital angular momentum. This equation can be easily integrated and the result is 
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p^ip — C -\-e/ pq, where C is a constant of integration that reduces to Lq in the absence of the 
perturbation. In fact, we can write C — Lq{1 + eA), which defines the constant A. We can 
now use the result of the integration of equation (5) to write equation (4) as a differential 
equation for = 1/p as a function of ip in the standard manner. The result is 

(Pu ^GM_ (3 + el) + 4eo cos (y? - Qq) 

The general solution of equation (8) is given by 

- _ 3 + ep _ CQip sin {ip - go) 



Loalil-el)^ Loag(l - eg)^ 

+ K cos {(p - go) + K' sin ((/? - go) , (9) 

where K and K' are integration constants. Let us note that equation (9) should reduce to 
the unperturbed ellipse given in equation (3) for e = 0. It follows from this requirement that 

K ^ /° , + efe , K'^ek' , (10) 
ao(l - e^) 

where k and k' are constants. It is possible to write equation (9) in the form 

1 + ecos [{l + S)ip-g] 

n T\ ' (11) 

a[l — e^) 



where (5 is a small dimensionless perturbation parameter defined by 

40; J cos i 



(12) 



Mc2(l-e2)3/2 ■ 

Here to — (GM/a^)^^'^ and the period corresponding to this "Keplerian" frequency is T = 
2t: /uj. Equation (11) contains three new constant orbital elements a, e and g defined by 

a = ao(l + ea), e = eo(l + e^), 5' = S'o + ^l^k' . (13) 

In a realistic situation, such as the motion of a satellite around the Earth, the gravitomagnetic 
perturbation caused by the rotation of the Earth cannot be turned off. Therefore, ao, eo and 
go should be thought of as auxiliary constants that are not directly measurable. On the other 
hand, by comparing the observed orbit with equation (11), one can determine the orbital 
parameters a, e and g. Using the definitions (13) in equation (11) and comparing the result 
with equation (9), we arrive at the following expressions for a, (5 and 77: 
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a 



1 



el 



2A + 



3 + 
Loao(l - el) 



+ 2eoaok , 



/3 = -(l-e^) , 
Co 

1 _ e2 



(14) 
(15) 
(16) 



1 + eg' 

where A is given by the relation C = Lo{l + eA) and k is defined in equation (10). 

The variation of (/? with time t is given by = C-u^ + eu^, where u{(p) is given by equation 
(11). Here C is different from L — ^GMa{l — e^); in fact, using equations (14) - (16) we 
find that 



C 



, 1 3 + e2 
L — -e 



(17) 



2 a(l-e2) ■ 

To solve equation (6) for Z, we employ the approach described in [21]: Let Z = pH{(p)] 
then, using equations (4) and (5) we can write equation (6) as 



Since H = 0(c~^) and FpH + F^p dH/d(p = 0(c~^), equation (18) implies that 



(18) 



dip^ 



sin ip + e sin ip cos {^p — g) + -e cos ip sin {ip — g) 



(19) 



where 5' is another small dimensionless perturbation parameter defined by 

, 4a; J sin i 





(20) 



Mc2(l-e2)3/2 ■ 

We assume that when observations begin at t = 0, (/? = (/9o and the instantaneous orbital 
plane coincides with the plane of the unperturbed orbit depicted in Fig. 1. That is, at 
</? = </?o we have H — and dH/dip — since both Z and Z vanish. With these initial 
conditions, the general solution of equation (19) is given by 

S' 

H — — [ sin + sin (ip — 2(^o) ~ 2((^ — (po) cos ip + e sin (ip + (po — g) 
+ e sin ((/? — (/7o) cos (2(/?o — g) — e cos (</? — (/?o) sin g — e sin (2(/? — g^) 
+ esin5r] . (21) 



This completes the construction of the orbit that is given hy X — p cos </? , Y — p sunp 
and Z — pH{(p). Here p — l/u{(p), where u{(p) is given by equation (11), and H{(p) is given 
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by equation (21). Prom the standpoint of observers in the basic {x, y, z) coordinate system, 
the orbit is characterized by the elements a,e,i,g,^} and (po, where Q is the longitude of the 
ascending node; indeed these six orbital elements can be determined in principle from the 
position and velocity of the sateUite ai t — 0. 

It is straightforward to show that the orbit thus constructed undergoes Lense-Thirring 
precession. To this end, one could start from the equation of motion (2) and study the 
average behavior of the Runge-Lenz vector and the orbital angular momentum vector as in 
[22]. Alternatively, one could study explicitly the orbit derived here. For instance let us note 
that for i — 0, the orbit is planar and the perigee occurs at (p — {l—S){g+2mr),n — 0, 1, 2, .... 
Thus the orbit precesses in the retrograde direction at the rate of 27rS/T — 4GJ/[c^a^(l — 
e^)^/^], as expected [19]. 

3. Azimuthal closure 

We assume that orbital observations are usually performed with respect to the fundamen- 
tal frame {x, y, z). The transformation between the {X, Y, Z) and (x, y, z) frames is given by 
the rotation matrix 



X 




cos VL 


—sin Vt cos i sin Vt sin i 




' X ' 




y 




sin f2 


cos Vt cos i —cos Vt sin i 




Y 


(22) 


z 







sin i cos i 




Z 





In the basic (a;, y, z) frame, let us introduce polar coordinates via x = r sin 9 cos 0, 
y = r sin ^ sin and z — r cos 9. Let us note that to 0(c~^), r = p along the orbit. We are 
interested in the time T that the particle would take to go from 0o at i = to 0o + 27r at 
t — T. The orbit is not closed in space; therefore, we seek the period for azimuthal closure 
with respect to the basic coordinate system that is used for observations. Since 

y sin cos ip + cos Q cos i sin ip — cos fl sin i Hip) , . 

tan0 = - = — ; — — — : — -— — . , , (23) 

X cos s I cos (p — smU cos ismp) -\- sm \ ls\n.i H [p) 

we evaluate tan0o from equation (23) aX p = p^ (corresponding to t = 0) and set it equal 
to tan (00 + 27r) evaluated from equation (23) at p — pr (corresponding to t — T). This 
simply implies, after some algebra, that for cos i ^ 0, 

sin {pq — pq) = [cos po H{pr) — cos pr H{po)] tan i . (24) 

Let us suppose that pr = po + 2% + A, where A is yet another small quantity; then, it 
follows from equation (24) that A = —n6' cos^ po tan i, since H{pq) = by assumption and 
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H{(pr) — — 7r5' cos (po according to equation (21). To ensure that |A| << 1, the inchnation 
angle i must not be near 7r/2 since tan (7r/2) = oo. Thus azimuthal closure in corresponds 
to </? : </?o — > V'o + 27r + A. 

It now remains to integrate ip — Cu^ + eu^ over the interval {ipo, (fir) in order to determine 
r, i.e. 

^"Zo Cu^ + eu^ ' ^^^^ 
where u{(p) is given by equation (11). To first order in A, equation (25) can be written as 

^ ^ dcp ^ A(l - e^/^ 

J(po Cv? + eu^ io\l + e cos{(po — g)]^ 

Moreover, equations (11) and (17) can be used to show that 

eip sin ((/? - g) 



uj[l-\-e cos (</? - g)Y 



1 + 25 



1 + e cos (</? — g') 



^1 1 — 2e cos (</? — 5') + 



(27) 



2 a;a3(l - e2)3/2 

Substituting this expression for the integrand in equation (26), we obtain T in terms of three 
basic integrals. To compute these, we recall that 

^ / s ^28) 



/co 1 + ecosC (l-e2)V2 ' J^^ (l + ecosC)^ {I - e^f/^ 
moreover, one can show via integration by parts that 

rCo+27r (^Q _|_ g^Q sin (^dC, TT TT 



ko (l + ecosC)=^ (1 + ecosCo)^ (1 - e2)3/2 ' 
The result can be written as 



(29) 



T GJ cos i j 3 4 — 2 cos^ (fQ tan^ i 1 

^ c^a^u; |~yr^+ [l + ecos((^o-5)]'J ' 

This expression reduces to our previous result [18] for a spherical orbit with e — 0. We 
emphasize that equation (30) properly contains only the gravitomagnetic correction to the 
period T to 0{c^^), i.e. T = T(l + Qgm)- The corresponding gravitoelectric perturbation, 
Qge, is given by equation (34) of the Appendix. Therefore, the final result may be written 
as 7I|- = T(l + Qgg ± Qgm), where the upper (lower) sign refers to a prograde (retrograde) 
orbit. 
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4. Discussion 

The prograde orbit for which we have derived the expression for T in equation (30) is 
such that the particle revolves around the source in the same sense as the proper rotation 
of the source (see Fig. 1). Imagine now the retrograde case such that at each point on the 
orbit the direction of the velocity vector is reversed. If at the same time the direction of J 
— and hence — is reversed as well, then + (v/c) x would be unchanged; therefore, 
the particle would be subject to the same force as before and relation (30) would hold in this 
case as well. It follows that if in Fig. 1 the orbit is assumed to be retrograde, then equation 
(30) would hold except that the overall sign of the perturbation term would be reversed. 
Thus the general formula for the clock effect is 



We expect that formula (31) would hold at 0(c~^) level as well if and 71 were replaced 
by the proper periods of spaceborne clocks on pro- and retro-grade orbits, respectively. 

It is interesting to note the dependence of the gravitomagnetic clock effect on the eccen- 
tricity of the orbit; this is illustrated in Figure 2 for an artificial satellite in orbit around the 
Earth. Though equation (31) is formally valid for any eccentricity, e cannot be too close to 
unity; otherwise, the small dimensionless perturbation parameters 5 and 5' given by equa- 
tions (12) and (20), respectively, would no longer be small and our perturbation analysis 
would break down. 

A remarkable feature of equation (31) is its topological character, i.e. its complete in- 
dependence from the semimajor axis of the orbit. Moreover, the gravitational constant G 
does not appear in equation (31). For the Earth, 47r J/(Mc^) ~ 10~^s; therefore, it may be 
possible to measure this effect in the near future [11-17]. 

The clock effect vanishes for a polar orbit. Equation (31) thus breaks down when the 
inclination angle i is sufficiently close to tx/2. This circumstance is related to the fact that 
for a spherical polar geodesic orbit, the orbit precesses with the Lense-Thirring frequency 
2GJ/ (c^a^) in the same sense as the rotation of the source. Therefore, the azimuthal closure 
period is in this case 27r{2GJ/c'^a^)~^. This period is extremely long compared to the orbital 
period; in fact, the corresponding ratio is ~ 10^° for near-Earth orbits. Thus in our pertur- 
bative approach i must be sufficiently different from 7r/2 such that |A| << 1. This fact is 
reflected in the requirement that in equation (30) the relative gravitomagnetic perturbation 
in the period be much less than unity. 




(31) 
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Finally, it is interesting to note that the clock effect (31) in general depends on the position 
of the particle on its orbit at t — 0. In fact, of the six orbital elements (a, e, i, g, Q, (fo), only 
a and Q are not represented in equation (31) due to the topological character of the effect 
and the assumed axial symmetry of the source, respectively. 
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Appendix 

The general solution of the gravitoelectromagnetic equation of motion to 0(c~^) can be 
obtained via a superposition of the perturbations due to £g and Bg. The latter have been 
the main subject of this paper. For the sake of completeness, we wish to indicate briefly in 
this Appendix the result of the perturbation analysis for £g. The work described in Sections 
2 and 3 can be repeated for £g as the perturbing acceleration. Indicating the corresponding 
quantities by a tilde, the perturbing accelerations are given by 



2 1 - Cq + 4eo sin^ v 
po ao(l - el) 



(32) 



F^p = — epo/Po = 0, where e = —AGMLq/c^. The analysis of Section 2 can be 

followed in much the same way with e — > e. The result is Z = and an orbit of the form 
(11) with S ^ S, where 

~r GM , , 

c^a(l — e^) 

and — Cu'^ + eu^ such that C is related to L as in equation (17) with e ^ e and 3 + ^ 
(5 + e^)/4. Equation (33) implies that the orbit undergoes Einstein precession by 2t\:\5\ 
radian per revolution. The period of the orbit T can be obtained from the integration of 
(f — Cv? + lu^ for ip : ipQ ^ ifQ + 2%, which agrees with (j) : <pQ —>■ (po + 2it except when 
i — 7r/2. The integration can be performed as in Section 3 and the result is T = T(l + Qge), 
where 



" 2c'a \ [l + ecos(^„-9)n 

For a circular orbit e = 0, we find O^e = SGM/ (2c^a). This has a simple physical explanation 
in terms of the isotropic radial coordinate employed in gravitoelectromagnetism. Imagine 
a test particle following a circular orbit in the Schwarzschild geometry. Suppose that the 
constant "radius" of the orbit is given by the Schwarzschild radial coordinate r. Let p be 
the corresponding isotropic radial coordinate, r = p[l + GM / {2c^ p)]^ . It is well known 
(see, e.g., [1]) that the orbital period in terms of r is given by the "Keplerian" formula 
27r(rVGM)i/2 = 27r(pVGM)i/2 [1 + GM/(2^p)f, which to 0{c-^) reduces to our result 
with p — a. 
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Figure captions 



Figure 1. Schematic plot of the unperturbed elliptical orbit. The orbital plane depicted 
here coincides with the plane of the osculating ellipse for the perturbed orbit at t = 0. 

Figure 2. AT = 71^ — 71 in seconds versus the orbital eccentricity e, < e < 0.95, for a 
satellite in orbit around the Earth. The common orbital parameters are (p^ = 45° and 
g = —195°. The solid line indicates an equatorial orbit {i = 0), while the broken line 
indicates an inclined orbit with i = 45°. 
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